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Abstract 

Thermal properties of the ordered phase of the spin 1/2 isotropic Heisenberg Antiferromag- 
net on a d-dimensional hypercubical lattice are studied within the fermionic representation 
when the constraint of single occupancy condition is taken into account by the method sug- 
gested by Popov and Fedotov. Using saddle point approximation in path integral approach we 
discuss not only the leading order but also the fluctuations around the saddle point at one-loop 
level. The influence of taking into account the single occupancy condition is discussed at all 
steps. 

PACS.75.10.-b General theory and models of magnetic ordering - 75.10.Jm Quantized spin sys- 
tems -75.50.Ee Antiferromagnetics 



1 Introduction 

The two-dimensional spin- 1/2 Heisenberg antiferromagnet (HAFM) on the square lattice has been 
extensively studied during the last few years. The motivation for this study stems from the 
discovery of high T c superconductivity in the ceramic compounds, where the competition between 
superconductivity and antiferromagnetic order has been observed experimentally M. 

Contrary to early suggestions there is nowdays strong evidence that the ground state of the 
fully isotropic quantum spin -1/2 HAFM on a two dimensional regular lattice is the Neel state (the 
classical ground state (the Neel state) is not disordered by quantum fluctuations.). 

This evidence is mainly based on numerical work 0] . Recently it has gained an additional sup- 
port by results obtained analytically with the help of various techniques e.g. large spin expansion, 
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field theory of the quantum nonlinear a -model || , effective Hamiltonian approach , a modified 
mean field approach (saddle point approximation) based on bosonic or fermionic representations 



The main problem in the technique based on these representations is to take into account the 
so-called single occupancy condition. 

The aim of this work is to study the thermal properties of the ordered (magnetic) phase of 
the spin 1/2 isotropic HAFM on a cZ-dimensional hypercubical lattice Q with periodic boundary 
conditions within the fermionic representation when the constraint of single occupancy condition 
is taken into account by the method suggested by Popov and Fedotov 0. We use saddle point 
approximation and discuss not only the leading order but also the fluctuations around the saddle 
point at the one-loop approximation level. We show that at zero temperatures one-loop corrections 
to the saddle point in our path integral description is equivalent to next-to-leading order in the 
linear spin wave theory. At all steps we discuss the influence of taking into account the single 
occupancy condition comparing the results of our calculations with those when this condition is 
disregarded. In particular we show that at finite temperatures taking into account the single 
occupancy condition considerably reduces the specific heat. 

For T ^ the two-dimensional spin system has no long range order (the Neel state is destroyed 
by thermal fluctuations) || and its state has to be treated as a paramagnetic one with strong 
antiferromagnetic correlations at finite distances. So our finite temperature results are relevant for 
the case when d > 3 and T < T N , where T N is the Neel temperature. 

In Sec. 2 we briefly review the fermionization procedure of spin operators by the method of 
Popov and Fedotov. 

In Sec. 3 we discuss the mean field result (the leading order of the saddle point approximation). 

In Sec. 4 we obtain the one-loop corrections (Gaussian fluctuations) to free energy and show 
that one can get the spin wave spectrum at zero temperture. We also find the specific heat and 
discuss the influence of the single occupancy condition on its temperature dependence. 

The last section is devoted to brief comments on our results. 

2 Fermionization by Popov and Fedotov' s method and bosonic 
path integrals for the partition function. 

The Hamiltonian of the isotropic HAFM reads 



the sum runs over ordered nearest neighbor sites of the d-dimensional finite regular lattice with M 
sites. For spin variables Si we assume periodic boundary conditions, J > . 

Many authors have proposed to use different representations of spin operators by Bose or Fermi 
operators. However, the fact that the dimensionality of the space in which these operators act is 
always greater than the dimensionality of the space of spin operators leads to the problem of the 
elimination of the superfluous states. Usually it is done by putting some constraints on the states. 

1 For simplicity we consider a simple hypercubical lattice though our approach may be also used for a non- bipartite 
lattice 



H of spin operators. 
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In the present paper we choose fermionic representation of spin operators 

Si=^4°-a^ i/9 , a,/3=l(1),2(D> (2) 

the summation with respect to repeated Greek indices is assumed, cr = (<j x , a v , a z ) are the Pauli 
matrices, and if we use this representation in the spin Hamiltonian H s we shall get the fermionic 
Hamiltonian: 

Hp = ^ {^^^{c^a^scjs) (3) 

<i,j> 

Ci a and c ia are fermionic annihilation and creation operators (at site i with spin projection a to 
the z axis), which obey canonical anticommutation relation 

Popov and Fedotov proved that the partition function of the model (1) 

Z = Tr s (e-^ s ) (4a) 

can also be written as 

Z = i M Tr F (.e- f>A *- i **). (4b) 
In these formulas Trg(Trp) is a trace in space where spin (fermionic) operators act, N = 

Y^tLi H a Cia = 2j=i ™* is t ne number operator. 

Let us briefly repeat their arguments. It is sufficient to consider only one site (we omit the site 
index). For spin i-, spin operators S a (a = x, y, z) act in two dimensional space. But the space 
where fermionic operators act is four dimensional; we have states 

|0,0), c t T |0,0) = |T,0), c t jO,0) = |0,|), c}c{|0,0) = | u>- 

States | T,0),|0, J.) can be identified with eigenstates of S z operator with spin up and spin 
down , we call them physical and denote \phys). Then states |0,0) and | f, J.) are superfluous or 
unphysical and their contribution should be excluded. 

The physical states span a two-dimensional physical subspace, characterized by the single oc- 
cupancy condition 

h\phys) = \phys). 

The direct product of the physical subspaces of all the sites form the sectors in which the 
Hamiltonians Hs and Hp coincide. 

In order to prove the basic formula Eq.(4b) we write 

Hp = Hpi + H' Fi , N = ni+Nl, 

where Hp^fii) is that part of the Hp(N) which contains the fermionic operators of the i-th site 
and H' Fi (Nl) is the remaining part. For the Hamiltonian of HAFM we have 

H Fi \unphys) l = 0. 

Therefore, the trace in Eq.(^) taken over unphysical states of the i-th site vanishes 
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since Tr, unp%s {(-i) ft< } = (-*) + H) 2 = 0. 

As a result , in the calculation of the trace all the unphysical states are eliminated, while 
the physical states Hf = Hs and N\phys) — M\phys). Therefore 



on 



Tr F (e 



) = (-i) M Tr phys (e-^n = ilVsfe-^) 



which proves Eg. fl4b|) . 

The evaluation of fcrmionic trace Tip requires only the standard technique because this trace 
is unrestricted. It can be represented as a path integral in terms of Grassmann fields r) and fj 



where 



Z = i M J T>fi v exp | - J 



Hf(v,V-, t ) = 



dr 



J 



(5) 



and 



4 E (Via {T)(T a pilif3 (t)) (?7 J7 (rja-ysViS (t) ) 

(i,3) 

J^S 4 (r).S,(r) 



(6) 



(7) 



0<T<P i,a 

Now let us do the Fourier transformation 



(8) 



qG-BZ m 



where q is the wave vector in reciprocal space (we can restrict it to the first Brillouin zone (BZ)), 
= is a Matsubara frequency for Bose field. 
Summation over ordered nearest neighbors can be written as 



E = ^E 



(9) 



(hj) 



since j is a nearest neighbor of i: Vj = + S, and 6 represents the displacement of z = 2c? nearest 
neighbors of each site . Then 



/ dTH F (v,V\T~) = J -^- Y E^ S ( q, ° m ) S ( -q ' -0 ™)' 

JO P „^D^ ™ 



(10) 



qEBZ m 
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where the so-called structure function 7q = -j J^S eZq ^ = 7-q = 3 (cos gi + cos (72 + ■ • • + cosqd). 
From Eq.(||) it follows that 

s(- q ,-n m ) = s*(q,n m ) (11) 

and if we write 

S(q O m ) = ReS(q, Q m ) + TmS(q, fi m ) (12) 

then 

/ drH F (fj, m T) = ^f - £ ^ 7q [(ReS(qfi m )) 2 + (ImS(q,n m )) 2 ]. (13) 

The standard way to decouple four fermion terms is to use Hubbard-Stratonovich represen- 
tations and introduce some auxiliary Bose fields. The decoupling scheme is not unique and the 
particular choice of the Bose fields depends which mean field solutions (ordered or disordered for 
our model) we are going to discuss. Of course, before one starts to use some approximation (usually 
saddle point approximation) all representations are equivalent and if we are able to calculate the 
path integrals exactly we shall get the same result. 

In the present paper we are considering the ordered phase (the disordered phase which is the 
most relevant for d < 2 will be discussed elsewhere M ) and the Hubbard-Stratonovich decoupling 
can be done with the help of an auxiliary vector field At(q, f2 m ) which plays the role of the 
staggered magnetization 



J lq,m 



V p 

and the path integration measure 



"'' M 7 q [M*(q, fi m )S(q, tlm) + S*(q, fi ro )Al(q, fi m )] 



, (14) 



^"=11 11 11 ~ ( 15 ) 



m q£_BZ a—x.y,z 

For Fermi fields Fourier transformations are 

ViajT) = } j ma{ l 'n)e~ WnT , (16a) 

n 

%a(T) = J2vMe iVnT , (16b) 

n 

where v n = 27r (^ +1 ) j s a Matsubara frequency for Fermi fields. Then 

S(q,n m ) = — Vi a { v ni)(T a i3ri l p{vn)e~ lq:r *8 nm . Unl ^ Vn . (17) 

i,n f ,n 

and for the partition function we get 
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Z = i 
where 

-a/3 



M 



^ c L- m i,j n' ,n J 



(19) 



and the measure for the path integration with respect to Grassmann variables now reads 

v ^n = II dViccivnjdriiaivn). (20) 



Integrating with respect to them we obtain 

Z = i M J V» M e- s °ffM =e -0F 



where 



and F is a free energy. 



Seff[M] = J2 \ M (^ Q rn)\ 2 ~ TvlnJC(M) 



(21) 
(22) 



q,m 



3 The leading order of the saddle point approximation 

In order to deal with the AFM solution we shall choose a frequency independent solution along 
the z-axis (tz = (ir, it, ir)) 



(23) 
(24) 



M(q,fi m ) = zy/dM0Jm5q t ir- 
The real parameter TYl is a staggered magnetization. Then 

(ic MF (M)r n ? n . !id = [s^ (-iv n p + i ^ + (-lydjpma^ 

and (— I) 1 = (— l)( ri ' lH i-( r i)d So the free energy in the mean field leading order takes a form 

F MF (m) = dJMm 2 - - ]Tln [1 + ef> Em ] mi, (25) 

where {E m } is the spectrum of the mean field Hamiltonian 



(26) 



3,a 
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u)jx = (-i) j djm, 
lu j2 = -(-i) j djm. 

The summation with respect to eigenvalues can be done easily with the result 

In [1 + e 0Em ] = M In ( ? cosh (d(3Jm)j . (27) 

So for the free energy in the leading order we get 

M M 

F MF {m) = dJMW? In (cosh (dfiJm)) In 2 (28) 

P P 

Minimization of F MF ( Tfl) yields the mean field staggered magnetization equation 

m = - tanh (d/3Jm). (29) 

Exactly the same result for magnetization one obtains in the mean field approach to the S = 1/2 
Heisenberg model with the Hamiltonian Eq.(Q) working in terms of spin variables. 

If the single occupancy condition is disregarded instead of Eqs.(pq) and Eqs.d2|) we get 



F MF (mo) = djMml - ™ in fcosh (^PJrrk) ) - ™ in 2 (30) 

and 



mo = itanh(^J772o ). (31) 



2 V 2 



4 One-loop corrections 

Now we write 



M(q, n m ) = zy/dMpjms^Tt + 8M(q, n m ) , (32) 

where (5A / t(q, VL m ) are fluctuations of the magnetization around the mean-field value (the leading 
order) 171 satisfying Eq.@. Then 

K',n-,iA 5M ) = (- iv "P + % \ ) + {-^djpmcjl^ 8 n , n 6ij + 5nf(v n , - v n ), (33) 
where 



qeBZ 

The partition function 



8nf{v n ,-y n )= \/-^7qe lq ' r ^Al(q,^ ~v n )cr a ^. (34) 
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Z = e-? FMF Jvn„e- s - a M BM \ (35) 

where 

S^ f [5M\=J2 \SM((q,a m ))\ 2 -Ti-\nJC(SM) , (36) 

m q£BZ 

X>Ma* = H 11 • (37) 

m qGBZ 

The superscript (2) means that only the terms of the second order with respect to 5M. are kept. 
Thus we take into account only so-called Gaussian fluctuations. 

Let us define a matrix G such that its matrix elements has a form 



7T 



n',n;2j 

This matrix is the one particle propagator evaluated at the saddle point. Then Eq.(|33|) written in 
the matrix form (with respect to spin index a and lattice site index i ) takes a form 

Kn> n = G~} n - 8n{v n , - v n ) (39) 

and 

TrlnK(SM) = TrlnG" 1 + Trln(l - G571) 

= Tr In G~ 1 - Tr [GS 71] - -Tr [GS TIG5 71] . 

The term which describes the Gaussian fluctuations in more explicit form reads 

Tr[G67lG5n] = ^ Sp[(m - w n /3)-Mn(fi m )(m - ii> n l3 + iQm^Sni-firn)], (40) 

where Vfl is a matrix in the space of i and a indices with the matrix element 

m°f = -(-l)°(-iydJ i 97?W tf *a/s, (41) 
and Sp is a trace in this space ( its element we denote as \ict)), 

(2) 

So in the one-loop approximation we have for S% f [8M] 

S tff = E ISMi^Qm^ + S^iSM] (42) 

and 

S$ f [8M] = ± ]T T aP (i,m){ia\Sn(n m )\iP){iP\67l(-n m )\ia), (43) 

i : ot,{3 : m 



s 



where 



= V - - (44) 

^ iP n P+ (-l) a {-l) l djpm iv n /3 - iflmP + (-l)P(-iydJj3m ' v 1 

The summation with respect to Matsubara frequencies can be easily done with the following 
result 



<$>(A,B;n m ) ^ 



1 



ii> n /3 - A iv n (3 - B - i£l m (3 

1 sinh[i(A-B)] 



in m f3-A + B i8mh[±(A + B)]+cosh[±{A-B)Y 

if A ^ B, and 

a a ^ x <■ 1 — isinhA 

A; fi m ) = -<5 m0 — -2—. (45) 

2 cosh A 

We need to know only the expression for the special choice: B = —A. In this case 

, , , ^ . 2AtanhA if2 m /3tanhA , . 

*^.-*.«~) = - l0 ^, + A A , - ia J ) , + tA , - <«> 

So we get (A = dJ/3m) 

m /■ \ r 1 - if-lV sinhA _ r w . . . 

-Tii(j,m) = -d m o 72—. = -Omo[K - i(-l) J p\, (47a) 

2 cosh A 

T22U, m) = -6 m0 [ K + i(-l) j p\, (47b) 

Tl2( ^ w) - - (sWT^ - (oi#T4^ s eM + C(m) (48a) 

and 

T21 (j,m) = C(m) - i(-iy((m). (48b) 
Taking into account the Eq.(f29|) for the mean field magnetization we get 

K = i(l-4m 2 ) . 

If the single occupancy condition is neglected we get instead (Aq = dJ(3Ulo ) 

*=;dhH (1 - 4TnS) 
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and 



e . , 2A)tanh^ TO /3tanh^ 



From Eqs. (p9| ) and (^Tj) it follows that £,{m)\A=A a — £o( m ) an d C( m )U=A = Co( m ) 
From Eq.(B4) we have 



(il\Sn(Cl m )\il) = ^ J-^ 7qe ^«5M z (q,O m ), (49a) 



i2|<Jn(fi m )|i2) = - ^ \/-^- 1(1 e^-6M z (q,n m ), (49b) 



il|5W(n ro )|t2) = 53 J-^W^lSM'faShJ-iSMifann)], (50a) 

qG-BZ 



q£BZ 



i2|jn(o m )|a) = 53 \/-^ 7 qe iqrt [^ a; (q,an) + ^ !/ (q I an)] , (50b) 



and <?ef < defined in Eg. ([43[) can be rewritten as 



where the longitudinal part 



5(2) _ 5(2) , 5(2) 

eff ~ L "T" J T ' 



4 2) = ^^H^ 1 ! 572 ^)!^^ 1 !^- ™)!^ 1 ) +7 1 2 2 M(*2|^n(r! m )|z2)(z2|«5n(-ar l )N2)} 

and the transverse part 

So using Eqs.fli^) and p9| ) we obtain 

§P=dJ0K 53 7q |<iM z (q,0)| 2 (51) 

and with the help of Eqs.(g§) and @ 

4 2) = 51 {7q[^ x ((q^ m ))-^^(q^m))][^"(-q,-an) + «^n-q,-« m )]CM 

+ 7q [£.M a; ((q, fi m )) + ttM^q, Q m ))}[5M x {-c L , -O ro ) - iW^-q, -O m )]e(m) 
-7q[^ a ((q, n m )) - iM^((q, fi m ))][*M a (-q - tt, -fi m ) + tfM»(-q - tt, -Q m )]C(m) 
^[^"((q, n m )) + i6M y ((d, n m ))][5M x (-q - tt, -Q m ) - iSM y (-q - tt, -fi m )]C(m)} . 
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So 



sfj } = dJpK l7ql{l^ 2 (q,o)| 2 -|5M 2 (q + 7r,o)| 2 } 

q£RBZ 

+ dJ PYl E W { [-I^ a (q ! m )| 2 + |^(q+7r,n m )| 2 (52) 

m qeRBZ 

\5M y (q, n m )\ 2 + \SM y (q + 7T, « m )| 2 ] CM 
+ 21m pM a (q, n m )$M»(-q - 7T, -fi m ) 
- 6M y (q,n m )6M x (-q-7V,-n m )]C(m) } . 

(2") 

The part of S^A which describes the transverse fluctuations for each q vector and for each 
value of m consists of four 2x2 blocks mixing the real and imaginary components of 5A4 X and 
SA4 y at q and q + 7r in pairs. The matrices corresponding to these blocks are given by (q € RBZ) 



1 - dj9J| 7q |e(m) ±d/3J| 7q |C(m) 
±d/3J| 7q |C(m) l + d/?J| 7q |£(m) 



(53) 



with the eigenvalues 



Mq,m) = l± , Mm|7ql (54) 

VOW" 1 ) 2 + (2dm) 2 

We see that A_ (q, m) vanishes at fi m = when q = . The corresponding eigenmodes are the 
Goldstone modes (spin waves) which appear due to the fact that AFM (ordered) phase is a phase 
with spontaneously broken symmetry. 

So from Eq.(p5|) we obtain up to some inessential constant the free energy including one-loop 
corrections 



F = F MF + - E ln (A + (q,m)A_(q,m)) + i £ ln(l- 7 2 ) 

P qeRBZ m>0 P qGRBZ 



l - Y ln(l - d 2 /3 W 7 2 ) 



2/3 



(55) 



qeRBZ 



where Yl'qeRBZ means that the point q = should be excluded. The summation with respect to 
Matsubara frequencies can be done with help of the formula 



EM 1 



m>0 



.4 2 



(fi m /3) 2 + B 2 



In 



B sinh 



VB 2 -A 2 



VB 2 - A 2 sinh f 



(56) 



and we get finally for the free energy 



F = F 



MF 



3 ^ 

H q€RBZ 



sinh 



(dmpjyf 



smh(dmpj) 



dmBJ 



-In 

3 n smh(dmBJ) 
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+ 



2/3 



1 



J2 ln[l-^ 2 J 2 (l-4m 2 ) 2 7 2 



(57) 



qERBZ L 



Now taking the limit of zero temperature — ► oo we get the energy of the ground state per 
site (there is no contribution to the limiting expression from the last term) 



which is exactly the ground state energy per site obtained in linear spin-wave approximation ]12[ 
for spin 1/2. 

It is easy to check that the same zero temperature result will be obtained when the single 
occupancy condition is disregarded. 

On the contrary, at finite temperatures taking into account the single occupancy condition gives 
different values of the different thermodynamical quantities, e.g. for the free energy this affects in 
changing the temperature dependence of the magnetization TTL and changing the longitudinal part. 
Figs.l and 2 show the difference of the results for the cases when single occupancy condition is 
taken into account (solid line) and when it is disregarded (dashed line). On Fig.l the temperature 
dependence of the internal energy and entropy and on Fig. 2 of the specific heat are given. From our 
numerical calculations we found that only in the interval < t < 0.13 the difference is negligebly 
small. In this interval the specific heat goes to zero as C v = at a with a = 3, as it should p2| . 

5 Conclusion 

We have studied the magnetic (ordered) phase of the isotropic spin 1/2 HAFM defined on the simple 
d-dimensional hypercubic lattice using fermionized spin operators and saddle point approximation. 
Single occupancy condition which is needed when spin operators are bosonized or fermionized is 
taken into account by the method of Popov and Fedotov. 

It is shown that inclusion of the one-loop corrections to the leading order of the saddle point 
approximation leads in the limit of zero temperature exactly to the same expression for the ground 
state energy which one obtains for the next-to-leading term in the linear spin wave theory, and 
this result does not depend if the single occupancy condition is disregarded or not. 

It is worthwile to mention that in the mean field theory of the spin 1/2 HAFM where Schwinger 
bosons are used [[| one obtains the same result of the linear spin wave theory at zero temperature 
already in the leading order and taking into account the single occupancy condition is crucial in 
this case. 

We demonstrated that in our approach taking into account the single occupancy condition 
changes finite temperature results considerably. 
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Figure captions 

Fig.l Internal energy per site E/M and entropy per site S/M versus dimensionlcss temperature 
t = (/3J) _1 for the 3-dimensional cubic lattice. Dashed and solid lines correspond to the cases 
when single occupancy condition is disregarded and taken into account respectively. 

Fig. 2 Temperature dependence of the specific heat per site C v /(KBt) for the 3-dimensional 
cubic lattice. Dashed and solid lines correspond to the cases when single occupancy condition is 
disregarded and taken into account respectively. 
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